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ABSTRACT 



An analytical derivation is pr esented for comput ing; mass-loss rates of 
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Cepheids by using the method of ICastor et al.l (119751 ) modified to include a 
term for momentum input from pulsation and shocks generated in the atmo- 
sphere. Using this derivation, mass-loss rates of Cepheids are determined as a 
function of stellar parameters. When applied to a set of known Cepheids, the 
calculated mass-loss rates range from 10~ 10 -10~ 7 M Q /yr, larger than if the winds 
were driven by radiation alone. Infrared excesses based on the predicted mass- 
loss rates are compared to observations from optical interferometry and IRAS, 
and predictions are made for Spitzer observations. The mass-loss rates are con- 
sistent with the observations, within the uncertainties of each. The rate of period 



change of Cepheids is discussed and shown to relate to mass loss, albeit the de- 
pendence is very weak. There is also a correlation between the large mass-loss 
rates and the Cepheids with slowest absolute rate of period change due to evolu- 
tion through the instability strip. The enhanced mass loss helps illuminate the 
issue of infrared excess and the mass discrepancy found in Cepheids. 

Subject headings: Classical Cepheids, Mass loss, Circumstellar Shells 



1. Introduction 



Cepheid variable stars are powerful tools for many aspects of astrophysics. Cepheids 
are excellent laboratories for stellar physics while also providing a period-luminosity relation 
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that is essential for extragalactic and cosmological studies. By observing both the intensity 
variation across the stellar disk as well as the variation of the angular diameter as function of 
pulsation phase, the emerging field of optical interferometry is extracting further information. 
This new information demands more sophisticated Cepheid models of the stellar interior 
and atmosphere, but also provides a more robust distance to individual Cepheids, helping 
to calibrate the period-luminosity relation that is more precise than any other calibration. 



Lane et al.l (120001 ) reported the first detection of the radius variation of a Cepheid, ( 



Geminorum, where the authors determine d the mean an g ular d iameter and incorporated 
radial velocity data to derive the distance. iKervella et al.l (120041 ) observed seven Cepheids, 
determining the variations of the angular diameter for four of them in the K-band with a 
precision better than 10%. In turn, these authors were able to calculate distances to these 
Cepheids with a statistical uncertainty of only a few percent, more accurate than HIPPAR- 
COS parallaxes. More recently, interferometric K-band observations of Polaris, S Cephei, 



and I Carinae provided evidence of circumstellar shells (IKervella et al.l 120061 ; iMerand et al. 



20061 ). These circumstellar shells cause an infrared excess of a few percent and increase the 
uncertainty of the angular diameter and the variations. The larger uncertainty, in turn, 
decreases the precision of the distance determination. The existence of these circumstellar 
shells pose a problem both for understanding the structure of Cepheids and for the calibration 
of a more precise period-luminosity relation. 



Merand et al.l (120071 ) continued to explore the phenomenon of circumstellar shells with 



interferometric observations of Y Oph and the non-pulsating yellow supergiant a Per. They 
found evidence of a circumstellar shell about Y Oph contributing approximately 5% of the 
total flux, larger than the excesses from previous observations. The yellow supergiant a Per, 
on the other hand, was found to have no flux excess. This result lends credence to the idea 
that the shells are related to radial pulsation. 

Interferometric observations appear to have raised a new set of questions regarding 
the evolution of Cepheids but the detection of circumstellar shells may be a manifestation 
of an older problem: the idea of mass loss in Cepheids. Throughout th e 1980's, mass 



loss was argued as bein g the cause of the mass discrepancy in Cepheids ( IWillsonl Il989 



Brunish fc Willsonl Il989l ). The Cepheid mass discrepancy (see ICoxi Il980l . for a detailed 
explanation) is the difference between predicted Cepheid masses estimated using evolutionary 
stellar models and tho se derived using pulsation models. The difference, at that time, was as 
much as 30% to 40%. IWillsonl (119891 ) argued there is significant mass loss related to shocks 
generated in the atmospheres of Cepheids. Furthermore the enhanced mass loss would slow 
the evolution of the Cepheids on the blue loop and increase the mass lost to winds even 
more. After some evolution, Cepheids could lose enough mass to reconcile the difference 
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between evolutionary and pulsational predictions of the mass of a Cepheid. Other theories 
were proposed to explain the discrepancy, such as convective core overshoot i n the progenitor 
models of Cepheid s, which could generate a steeper mass-luminosity relation (IKeller fc Wood 
20061 ; Keller! 120081 ) . Another theory proposes missi ng opacities in the evol utionary models. 



In the end, the introduction of the OPAL o pacities ([Rogers fe Iglesiaslll992l) had the greatest 
effect on reconciling the mass discrepancy; iMoskalik et al.l (119921 ) found the period estimates 
to be reduced for the same set of stellar parameters (L, M, T eS , X, Z) using the OPAL 
opacities compared to the Los Alamos opacities. This result seemed to end the theoretical 
study of mass loss in Cepheids as the current mass discre pancy is of the o rder of 10% for 
short period Cepheids, and decreasing with longer period (IBono et al.l 1200 ll ). 



While the theoretical study of mass loss in Cepheids is lim i ted, m uch more work has 



been presented on the observational side. iMcAlary &: Welch! (119861 ) used IRAS to ob- 



serve a sample of Cepheids to search for infrared excess. They found small infrared ex- 
cesses in twelve Cepheids and estima ted mass-loss rates of the order 10~ 9 — lO~ 8 M /yr. 
More IRAS observations (iDeasyl Il988l ) found upper limits of mass loss rates ranging from 
10~ 10 — lO~ 6 M /yr. These mass-loss rates were calculated by assuming the observed flux 
excess is caused by dust emission. Radio observations using the VLA set upper limits for 
mass-loss rates of order 1 0~ 9 — lO^M^/vr for the classical Cepheids FF Aql, rj Aql, SU 
Cas, 5 Cep and T Mon (IWelch fc Durid Il988l ) by assuming the emission is due to hot 



ioniz ed gas. With observations 



Rodrigues fc Bohm-Vitensd ( 



10m the IUE satellite of the Cepheid binary SU Mus- 



1992 f) determined an upper limit for the mass-loss rate 



cae, 

M < 7 x lO _lo M /yr. iBohm-Vitense fc Lovd (I1994J ) found a mass-loss rate for / Carinae of 
2 x lO _5 M /yr with an uncertainty of two orders of magnitude, again using IUE observa- 
tions. Collectively, these observations provide significant evidence for mass loss in Cepheids. 
Indirect evide nce of mass lo ss is the detection of circumstellar shells surrounding RS Pup 
(lHavlenlll972l ) and SU Cas ([Turner fc Evanall984l ) from optical observations, though both 
Cepheids are associated with reflection nebulae. 



The current state of the Cepheid mass discrepancy has been reviewed by IBono et al. 



(120061 ). where four plausible solutions are presented: extra convective mixing, rotation, stellar 
opacity and mass loss. If one increases the amount of convective core overshoot in main 
sequence evolutionary models of Cepheid progenitors then more hydrogen mixes into the 
core and creates a more massive helium core. This causes a more luminous Cepheid at fixed 
mass. Increasing the rotation in Cepheid progenitors has the same effect as extra convective 
mixing while increasing the radiative stellar opacity would affect the driving of pulsation, 
causing a longer period for the same mass or conversely a lower mass for the same period of 
pulsation. Because new stellar opacities went far in resolving the mass discrepancy before, it 
is possible the discovery of new radiative opacities at temperatures 10 5 K to 10 6 K could still 



-4- 



make a difference. Mass loss is also a possibility, and the connection between mass loss and 
pulsation has been explored recently for other type s of stars, suc h as lumino us blue variables 
(|Guzik et al.ll2005f ). asymptotic giant branch stars (jHofnerlll999l ) and Miras (IHofner &: Dorfi 
1997h. 



If circu mstellar shells surroun ding Cepheids are generated by mass loss then the ob- 
servations of iMerand et al.l (120071 ) imply mass loss is driven by pulsation, however, this 
does not preclude other driving mechanisms. The most common mechanisms for driving 
a stellar wind are: dust dr iving, radiative line-driving, coronal winds, and magnetic fields 
( ILamers &: CassineUilll999l ). Dust driving is important in stars where dust can form in the 
atmosphere. This leads to an increase of the opacity of the star and radiation pressure ac- 
celerates the dust away from the star. If the dust is ionized then the particles have a larger 
effective collisional cross-section and causes more material to be ejected. Since dust must 
form in the atmosphere for this to be important, the star's effective temperature must be 
cooler than 2000 K; Cepheids have effective temperatures of order 4000 K - 6500 K. Dust 
is not a likely factor though it can form some distance from the star where the temperature 
of the gas is much less when the gas is optically thin. 

Radiative line-driving has been explored thoroughly for most stars and is considered to 
be most important for hotter stars where carbon, nitrogen and oxygen atoms are ionized in 
the atmosphere causing a larger effective collisional cross-section. The ions are accelerated 
by radiation and interact with other atoms in the atmosphere, in turn accelerating them 
outwards. For stars with an effective temperature of order 6000 K, neon, iron, h ydrogen 



and h elium may play the role that carbon, nitrogen and oxygen play in hotter stars ((Abbott 



19821 ). Therefore radiative line-driving must be considered in any analysis of mass loss in 
Cepheids. 

T o have a coronal wind, a star must have a corona whe r e the temperature is very 
large. ISasselov fc Lester! Jl994f ) and ISchmidt fc Parsons! Jl982l . ll984aB ) found only weak 
chromospheres for Cepheids implying there are no hot coronae; coronal driving is an unlikely 
source. 



The presence of magnetic fields is uncertain and c ontroversial in Cep heids. iPlachinda 



( 120001 ) detected a 100 G magnetic field in 77 Aquilae but IWade et al.l (120021 ) could not verify 
that result, inferring an upper limit a 10 G field. However, Polaris has recently been observed 
to be a soft X-ray source, though this may be due to its binary companion; far UV obser- 
yations of Polaris imply the possible existence of warm winds, shocks or magnetic activity 
( lEngle et al.ll2006l ). Magnetic fields seem an unlikely cause of winds in Cepheids. Thus one 
needs to be concerned with only one of those four possibilities: radiative line driving. 
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The purpose of this study is to explore the role radial pulsati on plays iu driving mass 
l oss in classical Cepheids when coupled with radiative line-driving. ICastor. Abbott. &: Klein 
( 119751 ). hereafter CAK, devised a method to solve for the radiation-driven wind structure 
of a star. Here that method is modified to include a simple function describing additional 
momentum input into the wind. The second section describes the derivation of the CAK 
momentum equation for a pulsating star, including pulsation and shock physics. Section 
three explores the mass-loss rates for observ e d Cepheids, and Sect i on 4 c omp ares predicted 



flux excesses with those determined by lDeasyl (119881 ) . iKervella et al.l (120061 ) and lMerand et al. 



( 120061 . 120071 ). The fifth section will consider the dependence of the rate of period change, 
P/P, on mass loss. 



2. Pulsation Enhanced CAK Method 



Castor et al.l (119751 ) proposed an isothermal wind model to describe the mass-loss rates 



and wind structures of stars. The method is a powerful, yet simple, analytical tool for 
understanding the effects of radiative line driving because it requires knowledge of only the 
mass, radius and luminosity of the star. The authors assume the forces due to the radiative 
lines and the pressure gradient are functions of the local velocity gradient; this results in 
having the conservation of momentum written in a form that can be solved for the velocity as 
a function of distance from the stellar surface. The CAK wind model is the standard tool for 



F and G supergiants (ILamers fc Cassinelli 



found in ILamers fc Cassinellil (119991 ). 



understanding winds in O, B, and A stars (Owockill2005l ; IWatanabe et al.ll2006l ) and possibly 



19991 ). A review of the analytic derivation can be 



2.1. Derivation for Pulsation Driven Winds 



One can derive the equation s governing the mass- l oss ra te for a radially pulsating star by 
following the derivation given by ILamers fc Cassinellil (119991 ). who start with the momentum 
equation for the wind including the effects of continuum and radiative line driving, gravity 
and gas pressure, 

GM* I dp 

+ --T+9e+gL- (1) 

p dr 



dv 



dr r 2 

The functions g e and gi describe the radiation force per unit mass due to electron scattering 
and line opacity respectively. The expression for g e is 



9e 



(T e L* 

47rr 2 c 



(2) 
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where T e = <j e L*/ (AncGM*), and a e = <TriV(e) is the electron scattering opacity as a function 
of the Thomson scattering opacity and the number of electrons. The electron scattering term 
is clearly dependent on the composition and temperature of the star, but for this analysis 
is assumed not to vary significantly. The acceleration due to lines is found by using the 
Sobolov approximation; the final form is 



9l 



Cr 



r 2 v 



dv 
dr 



where 



C 



Aixc Z 







CT e VthM 
47T 



10 



-ii 



W 



(3) 



(4) 



The function C is treated as a constant, although this is not strictly true. The parameters 
a, k and 5 are force multiplier parameters calculated using model atmospheres to measure 
the radiative acceleration due to lines. The values of a and 5 tend to be of the o rder 0.5 
and 0.1 respectively, while k can vary significant ly: for stars T e g < 10 4 K k « 0.1 ( lAbbott 



19821 ; IShimada et al.lll994j ; iPauldrach et al.lll986l ). The ratio of the electron density with the 
geometric dilution factor W(r) = (1 — \/l — R*/r) /2 can vary for all stars implying C is not 
constant, but the ratio has the exponent 5 that is small. Therefore the term to the power 5 
is ignored by assuming it is approximately unity. It should be noted, however, in Cepheids 
the ionization rate of hydrogen and helium is dependent on the phase of pulsation. The 
number density of electrons could thus vary as a function of phase and play a small role. If 
the number density of electrons is increased then the mass-loss rate that is calculated may 
be larger as a result, but not by more than an order of magnitude because the term T e is 
much less than unity. Also C has an explicit and implicit metallicity dependence found from 
calculations on the contributions to radiative line driving for variou s ions as a function of 
gravity and effective temperature (lAbbottlll982l ; IShimada et al.lll994l ). 



The pressure term is rewritten in terms of the velocity by assuming the wind is isother- 



mal 



1 dp 
p dr 



2a 2 dv 2a 2 
r 



(5) 



v dr 

where a is the isothermal sound speed. The resulting total form for the conservation of 
momentum is 

a 2 dv 2a 2 
v dr r 



dv 



GMJi-r r 



V 



dr 



+ — ^ + — + Cr~ 2 [ r 2 v— I 



dv\ 

' — 
dr ) 



(6) 



Equation [6] does not include the effect of pulsation. Radial pulsation injects momentum into 
the wind due to the acceleration of the outer layers of the atmosphere of the star. Radial 
pulsation also generates shocks in the interior that propagate to the surface, depositing energy 
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into the wind. This means there are two additional sources of acceleration in the atmosphere 
of the Cepheid. For this derivation, it is assumed the accelerations due to pulsation and 
shocks are determined at the surface of the star so that they are only a function of time 
and global parameters. In the wind, these additional acceleration terms are modified by 
a dissipation factor, (r/R*)~ u where v > 0. Therefore the sum of the acceleration due to 
pulsation and the acceleration due to shocks are written as ((r/R if )~ u . The function ( is 
the acceleration at the surface of the star due to the sum of the pulsation and shocks that is 
written in terms of the global quantities of the star. The function is derived in detail in the 
next section. The full equation for the conservation of momentum including pulsation is 

dv GMJl — T e ) a 2 dv 2a 2 „ 2 / „ dv\ a , f r \~ u 
rv— = *\ - + — — + + Cr- 2 [r 2 v—\ +([ — ) ■ (7) 



dr r 2 v dr r \ dr J \R* 

Rearranging terms and simplifying, Equation [7] becomes 

c d^i) " - (> - $) ^ - GM ^ ~ r <> + ^ + (i) M = »■ < 8 > 

Equation[8jrepresents the quasi-static approximation of the momentum equation for the wind 
of a pulsating star. The differential equation can be solved in a series of snapshots where the 
velocity calculated is based on the conditions of the star at the given time, however this does 
not change the mass-loss rate at that instant. The mass-loss rate depends on the velocity 
and the density of the wind as it is ejected, but the wind does not feed back onto the layers 
that generate the wind. Therefore the time averaged mass-loss rate can be determined by 
this quasi-static approximation but the velocity of the wind will be unphysical near the 
surface of the star. Hence only the change of the luminosity and radius and other quantities 
dependent on them are considered in the problem. 

Because the implicit effects of time are ignored, the differential equation can be solved 



in the same manner as the non-pulsating form. From Figure 4 of ICassinellil ( 119791 ). it is 
clear there is only one solution to the non-pulsating differential equation, the solution where 
the velocity as a function of distance from the star satisfies both regularity and singularity 
conditions at some critical point, r c . It must be noted the quantity r c is a mathematical 
construct for solving the differential equation and is not necessarily a physical quantity. If 
F — represents the momentum equation, Equation [HI then 

'dF\ 

— — = Singularity Condition (9) 

dv'J c 

( dF \ f dF \ 
— + ( v r -?r— ) =0 Regularity Condition, (10) 



\dr J \ dv 

where / denotes the derivative with respect to r. While it may seem that using the velocity 
structure to determine a valid solution is contradictory with previous statements, this would 
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be a reasonable approximation if the star were non-pulsating. The amount of mass ejected at 
time t is not be affected by the amount of mass ejected an instant later. At that instant when 
the mass is ejected it has a velocity structure for a steady flow that provides a satisfactory 
solution to the differential equation, though this ignores the effect of clumping. The singular 
and regularity conditions become 

[l-^)rlv c v' c = aC{rlv c v' c )\ (11) 



2aC 



+ - 



R* 



l-v 



aC(r 2 c v c v' c y 



2^r 2 c v c v' c 

V c 



(12) 



respectively, where one should recall d(/dr = as ( is defined at the surface of the star. 

One may now use the momentum equation and the singularity and regularity conditions 
to determine the velocity and velocity gradient at the critical point by first substituting the 
momentum equation into the singularity condition, and eliminating the quantity C, 



a 



a 



i — — x ' 



GM*(l-T e )-2a 2 r c -(Rl(^- 



(2-v) 



(13) 



Combining the singularity condition with the regularity condition will determine the velocity 
gradient at the critical point and using that result with Equation [13] one finds the velocity 
at the critical point 



1 + 



2^2 



(l-v) 



1/2 



(14) 



2 2 , 

v r = a + 



1 + 
GMJl-T 



2a 2 V R* 



l-v 



1 -1/2 



a 



1 — a 



x 



2a< -(R*(-^- 



l-v 



(15) 
(16) 



These two expressions give the velocity and velocity gradient at the critical point and form a 
starting point for solving the momentum equation. The expression v' c is combined with the 
singularity condition to determine the radiative driving constant C in terms of the radius 
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and velocity at the critical point and various terms from the momentum equation 



C 



1 — a 



1 CBj^-v) f T 



2a 2 



R* 



GM*(l-r e )-2a 2 r c -(Rl( 1 j r 



(l-v) 



2-v 



-a/2 



/p(X^>20L 

c u c 



(17) 



Equating the two expressions for C, Equation HI and [T7] and using the function for the critical 
velocity, Equation [T5J, yields the relation for the mass-loss rate 



M 



a e L*k Z 

47TC Z Q 



1 



[I- a) 



Air 

& e Vth 



GM*(l-r e ) -2a 2 r c -(Rl 



2 1 ' c 



R* 



2-v 



-1/a 



x < a r r 



1 + 



CR*{2-v) ( r 



2a 2 



R* 



(i-O' 



1/2 



+ 



a 



1 — a 



GM*{l-T e )-2a 2 r c -CRl[jj- 



2-v 



18) 



It is now shown the mass-loss rate of a Ceph eid can be calculated by k nowing the global 
parameters for the star and the critical point. lLamers fc Cassinellil (119991 ) argue that if one 
assumes a value for r c then the velocity structure can be solved and the density structure can 
be derived from the velocity. The continuum optical depth of the wind is calculated from 
the density. This lends itself to a condition for the solution of the differential equation. The 
stellar surface is generally defined where the mean optical depth equals 2/3, therefore the 
optical depth of the wind, given by the integration of the density and opacity from infinity 
to the stellar surface, must be 2/3. The inclusion of pulsation and shocks, however, means 
it is not so clear this criterion needs to be satisfied for a wind to be driven. Because there 
is no ideal criterion to use that will give a clear solution, the optical depth r = 2/3 will be 
used. The mass-loss rate is tested and shown to verify that it varies only weakly for different 
values of the critical point; hence it is reasonable to use the integral of the optical depth as 
a criterion for solving the momentum equation. Therefore at some phase of pulsation, one 
can calculate the instantaneous mass loss rate by adopting a value for r c and solving the 
wind structure and comparing the velocity predicted at the surface with the velocity due to 
pulsation; a solution is reached when the two values match. 
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2.2. Defining the Acceleration Due to Pulsation and Schocks 

In deriving the solution for the enhanced mass loss due to pulsation and shocks, it was 
assumed the function (, which represents the sum of the acceleration from both pulsation 
and shocks, depends only on the global parameters that describe the Cepheid, such as 
effective temperature, radius, amplitude of radius variation and period of pulsation. In this 
subsection, the formulation for the acceleration will be defined and justified. 

For simplicity it is advantageous to consider a one-zone model for the radial pulsation 
of the Cepheid with period P. In that case the radius of the Cepheid as a function of time is 
written as AR(t) = — AR cos(W), where AR is the amplitude of pulsation. The amplitudes 
of the neares t Cepheids are readily availa ble from radial velocity studies, and a list has been 



compiled by iMoskalik fe Gorynyal (120051 ) . The minus sign in the equation defines a phase 
of zero to correspond to the minimum radius and u = 2tc/P is the angular frequency of 
pulsation. The acceleration of the surface of the Cepheid due to pulsation is thus 

a puls = d 2 [AR(t)]/dt 2 = Lj 2 ARcos{ut). (19) 

The luminosity varies as a function of time as well. Since there is a well-known phase 
lag, such that the maximum luminosity occurs roughly a quarter period before the maxi- 
mum radius, the quasi-adiabatic approximation is used to describe the luminosity such that 
AL(t) = AL sm(ut). The electron scattering opacity can be written via Kramer's law to 
be a e oc pT~ 3,5 , implying it too varies due to pulsation because both the temperature and 
density vary as a function of time. Because the temperature goes as T c g- oc [L*(t) / (R^t) 2 ] 1 /* , 
the variation of the temperature is 

ATeff(t) _ 1 AL*(t) 1 ARJJ) 
T e g 4 2 R* 

In addition, the linear perturbation of the density in the one zone model is 

M = _ 3 H^, (21) 
P R* 

meaning the variation of the electron scattering opacity is 

Aa P 5 AR 7 AL 

(22) 



where variables denoted with a bar are the mean values of the variable over a pulsation 
period. Also, because the effective temperature will vary as a function of phase then so will 
the isothermal sound speed, which goes as a oc T 1 / 2 . This defines all of the quantities that 
are assumed to vary due to pulsation. 
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The role of shocks is not as easily quantified as the re is no s imple met hod to describe 
the ac celeration of gas by spherically symmetric shocks. IWillsonl (119761 ) and IWillson &: Hill 
(119791 ) develop a model for a shock in the at mosphere of long per iod v ariables where the 
shock is periodic with the pulsation. However iMathias et al.l (120061 ) and ISasselov fc Lester 
(jl990l ). for example, find evidence for multiple shocks per period. Also, while the mechanism 
for generating shocks is understood to be related to both the opacity mechanism and the 
7-mechanism, the strength of shocks is uncertain. An analysis of shock behavior in 6 Cephei 
was conducted by iFokin et al.l (119961 ). who modeled the atmosphere of S Cephei and traced 
the difference between the post-shock velocity of the gas and the pre-shock velocity of the 
gas in the frame of reference of the traveling shock as a function of the phase of pulsation. 
Our work uses that information to estimate the acceleration of the gas in Cepheids due 
to shocks. Here it is assumed that the velocity of the shock itself goes to zero at the 
surface of the Cepheid. This assumption is justified based on evi dence that Cepheids have 
weak chromospheres generated by shocks (ISasselov fc Lesterlll994j ). Therefore the difference 
between the velocity of the pre- and post-shocked gas at the stellar surface in the frame of 



the shock is the same as in the rest frame. It is also assumed that the velocities in Fokin et al. 



( 119961 ) are valid in a one zone model, i.e. these velocities are at the surface of the Cepheid. 



From iMatzner fc McKed (119991 ) and iKlimishin fc Gnatyki (ll98ll ) , one can assume the 
pressure of the post-shocked gas is proportional to the mean energy density 



u = — 



E 

m ' 



(23) 



where u, E, m represent the velocity, mean energy density and the mass of the post shock 
gas, respectively. The mean energy density is roughly constant, while the mass is not. For 
a central explosion, such as a supernova, one can write the mass as m « pr 3 , where r is 
the scale of the post shock gas from the shock. However, shocks in Cepheids are generated 
relatively near the surface, therefore we can write the mass asm^ pRlAr, where Ar is the 
thickness of the shock front. Thus the velocity relative to the sound speed inside the star, 
c s , scales as 

Co 



(24) 

where Q is a constant, and Ar is the mean free path times the pre-shock pressure divided 
by the change of pressure (jZel'Dovich fc Raizerlll967l ). 



Ar = / 



V 

Ap' 



(25) 



The pre-shock pressure is c 2 s p at the surface and the change of pressure is the sound speed 
at the radius where the shock is formed times the change of density. The shock is formed 
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at the partial ionization zone, which has the same temperature for all Cepheids. The mean 
free path is 1/np, and the opacity of the shock is also constant. Therefore the thickness of 
the shock front is 

Ar = — eS = — * R * = —— (26) 

K ApT ionization «ApLV4( r ) r -i/2 K A Pi? y 2 ' 1 ] 

where the luminosity is the same at both r and R*, and the radius where the shock is formed 
is roughly constant for all Cepheids. Thus the speed of the gas due to a shock, given by 
Equation [2U is proportional to 

u oc c s (pRTf 2 . (27) 

Therefore the ratio of the velocity of gas for any Cepheid relative to the prototypical Cepheid, 
S Cephei, is 

u c s p 



7/r -l/2 D 3/4' 

us c S)& p 6 > R/ 



(28) 

where 6 refers to 5 Cep. The sound speed is proportional to the square root of the effec- 
tive temperature which, in turn, is proportional to the luminosity and the radius, T e fj oc 
L 1 / 4 /?^ 1 / 2 . The relative shock speed is now u/us = (L ^k 1 / 8 pR^J 2 ) / (L^ 8 'pgRg ) ■ The radius 
can be expressed in terms of the mean density and the mass, R* oc p~ 1 / 3 M 1 / 3 while the 
mean den sity is a funct ion of the period of pulsation p = (Q/P*) 2 where Q is the pulsation 
constant. iFernid (119671 ) found that pulsation constant is not strictly constant, and actually 



u s [lsJ \Ps) \Ms) 



varies as Q oc P 1/s . Equation [28] is rewritten in terms of the luminosity, period and mass, 

(29) 

This defines the velocity of the gas due to shocks at the surface of a Cepheid. 
The acceleration of the gas is a s h oc k = P~ 1 du/d(j). In terms of Equation | 

/M v 7nv 7/2 V^V /12 



ashock = —du s [ — j (-] [ w \ , (30) 



where the quantities describing 5 Cep are taken from the model computed by iFokin et al. 



( 119961 ). The values of du/dcj) are listed in Tabled], where zero phase corresponds to minimum 



radius. 

Having defined the acceleration of the gas due to shocks and pulsation, the function £ 



is 



- 'M 1/8 fp.y V2i (ma s/12 



C = A^oosM) + ^ Ui ^J (j^J , (31) 

where ( is dependent only on global parameters: the amplitude of radius variation, luminos- 
ity, period, and mass. 
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Range of Phase 


du/d0 (km/s/phase) 


- 0.08 


187 


0.1 - 0.12 


208 


0.13 - 0.22 


120 


0.65 - 0.68 


75 


0.73 - 0.76 


200 


0.79 - 0.83 


240 


0.93 - 1 


187 



Table 1: The change of speed of shocked gas at the surface of 5 Cephei with respect to phase 
of pulsation at different phases. At any phase not within t h e give n ranges the change of 
speed is zero. The data are taken from Figure 4, iFokin et al.l (119961 ) . 



2.3. Solution Space of the Wind Equation 

It is important to explore the momentum equation's parameter space before applying 
it to real stars. There may be situations where the combination of parameters do not give a 
physical solution, which is true for the momentum equation for static stars. There are two 
regions in the velocity-distance parameter space where there are no solutions: one region is 
defined where the velocity is subsonic when GM*(1 — T e ) — 2a 2 r < 0, the other is when the 
velocity is super sonic too close to the surface of the star where GM„(1 — T e ) — 2a 2 r » 0. 
Cassinellil ( 119791 ) reviewed the solution space of a radiative wind. The first case is e xplored 



by analyzing how the addition of pulsation affects the Parker point, r p ( lParkerlll958l ). which 



is the distance from the star where the effective escape velocity equals the isothermal sound 
speed. In the pulsating case the Parker point is defined as 



GM*(1-T e ) -2a 2 r p -(Rl 



2-v 



0. 



(32) 



Second, the boundary where the wind becomes supersonic too close to the surface will be 
explored by analyzing the parameter space of the critical velocity and the derivative of the 
critical velocity. 

The Parker point is an important quantity in the radiative driving momentum equation 
as the wind solution must be supersonic at a distance less than 



r < TV 



GM»(1 -r e 
2a 2 



(33) 



otherwise the velocity gradient dv/dr < and the v(r) would go to zero. The inclusion of 
pulsation adds an extra complication since the Parker point is now defined as the solution 



-14- 



to a polynomial equation. In the momentum equation given by Equation [HJ it is necessary 
that 

F(r) = GM»(1 - T e ) - 2a 2 r - QR v y~ v > (34) 

when a 2 /v 2 > 1. Because this function is non-linear, the Parker point is not a complete 
definition of the instability, especially since it is mathematically plausible to have a solution 
where r p is less than the stellar radius i?*. The true Parker point can be defined by using 
F{r p ) = and dF(r p )/dr < 0. The second definition ensures F(r) is decreasing before the 
Parker point. This is seen if different values of v are assumed; consider v = 0,2, 3, which 
are the quadratic and linear cases for Equation [32] for the Parker point. In the linear case of 
v — 2 the Parker point is 

In this case the Parker point is modified from the traditional definition to one that oscillates 
about GM»(1 — T e )/2a 2 . This case does not suffer from multiple solutions, but the cases 
^ = 0,3 have two solutions each, where 



and 



for„ = 0, rp= 2a^^ + 4CGM,(l-r,) ^ 

_ gm.(i - r e ) y[GM4i - r e )p - mm m) 

tor v 3, r p - — ± — [67) 

respectively. These solutions highlight another potential pitfall, the term in the square root 
may be imaginary depending on the phase of pulsation. In the case of v = 0, if the phase 
0.25 < <p < 0.75 then it is possible for Q < when the shock acceleration is zero and the 
acceleration due to pulsation is less than zero. For v = 3 the Parker point can be imaginary 
when the function ( is large as is the case for short period Cepheids or hotter evolved 
Cepheids. If the Parker radius is imaginary then there may not exist a wind solution at that 
time as the gravity is less than the outward acceleration implying an instability in the star. 

Applying the second criterion for a proper solution for the Parker point, one finds the 
derivative of Equation | 



dF{r p ) 



2a 2 - (2 - v)(R* ( If ) . (38) 



X-v 



dr y \R 

Testing this for the case where v = 2 the derivative is always less than zero. When v = 0, 
the Parker point is given by Equation [36] and substituting this result to solve for dF(r p )/dr 

= ±^ + (GM»(l-r e ) (39) 
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implying only the negative square root in the solution for the Parker point as given in [37J is 
the true Parker point. From Equation [371 it is clear the Parker point given by the positive 
square root is less than for most values of £. 

The second region of interest is based on the wind becoming supersonic too close to 
the surface of the star, where the boundary is defined as the curve that is tangential to the 
critical velocity. Therefore Equation [15] can be used to analyze any instabilities. There is 
the obvious possible instability given by the term {l + [(( y 2-v)R Jf /2a 2 ]{r c /R Jf ) 1 - u }- 1 / 2 . This 
term can be imaginary depending on the chosen value of v and the phase of pulsation. If 
this is imaginary, then both the critical velocity and the derivative of the critical velocity 
is complex and the will be no real wind solution. It should be noted that this instability 
disappears in the case of v = 2, and it is possible to compute a mass-loss rate at any phase 
of pulsation. 

The conclusion is that there exists time-dependent perturbations to the Parker point 
and the critical velocities causing no real wind solutions. This is different from the non- 
pulsating CAK method which does not contain any explicit pitfalls. 



2.4. The Power— Law Dependence of the Pulsation and Shock Acceleration 

It is important to understand the role the exponent v plays in the pulsation term of 
the wind equation. The power law is chosen to represent the dissipation of energy imparted 
by pulsation as particles in the wind are accelerated outwards. The use of a power law is 
an approximation of the detailed physics of the dissipation process. Without having a firm 
physical determination of v, it is necessary to determine how sensitive the total mass-loss 
rate is to an adopted value of v. One way to test the effect of v is to derive the change 
of instantaneous mass-loss rate with respect to v. This will be done by calculating an 
approximation of the instantaneous mass-loss rate and solving it for various values of v. 

The sensitivity of the instantaneous mass-loss rate with respect to v can be measured by 
the quantity d\a.M /dv. The mass-loss rate is given by Equation [18], which can be rewritten 
in the form M = K\A l l ^ a (A 2 + K 2 A\); the term K\ and K 2 = a/(l — a) are constant 
multipliers, and A% and A 2 are the respective functions of v. The derivative of the mass-loss 
rate is 

dlnM 1 1 dA x dA 2 l ' dv + K 2 dA 1 l ' dv , . 

= H - — (40) 

dv a Ai dv A 2 + K 2 A X v ' 

where 

/ \ 2 ~ v 

A l = GM*(l-r e )-2a 2 r c -(R 2 J^) , (41) 



R* 



-16- 



A 2 = a 2 r r 



CR*{2-u) ( r. 



ul V2 



2a 2 \R 

This can be considered one term at a time; starting with A\ 

dA 1 



. + 



(42) 



dv 



and 



1 1 dA l ln(r c /R*) 



a A\ dv a 
The derivative of the second term is 



-1 + 



GM*{l-T e )/r c -2a 2 



(R* (r c /R* 



,1-1/ 



-1 



(43) 



(44) 



dA 2 
dv 



a 

y 



C(2 - v)R, ( r, 



2a 2 



R* 



(i-i/) 



-1/2 



(R* ( r. 



2a 2 \R 



(l-u) r 



1 + (2 - z/) In 



R* 



• (45) 



Starting with the result for dAi/dv, it can be seen that there is a large change with respect 
to v if the pulsation term is of order the effective potential, where the enthalpy term given by 
the isothermal sound speed is relatively insignificant. In the limit the pulsation term is equal 
to the effective potential then the ratio approaches one and the term in square brackets goes 
to zero. Therefore dA\jdv — > — oo and the derivative of the mass-loss rate will approach 
infinity. The terms (dA 2 /dv) /(A 2 + K 2 Ai) and K 2 (dAi/dv)/(A 2 + K 2 Ai) can be shown 
to have a similar behavior. Therefore v can have a significant effect on the mass-loss rate 
depending on the combination of parameters. 

While the choice of v affects the mass-loss rate, it also affects the predicted density 
structure that is used to calculate the continuum optical depth as a criterion for a solution 
to the momentum equation. Consider the mass loss for a pulsating star from two different 
dissipation laws given by v\ and v 2 where v 2 > v\. For the power-law v\ there is a solution 
given by a corresponding critical point r c l , which satisfies the Singularity and Regularity 
condition, Equations M and [TD] respectively, as well as predicting a density structure such 
that the continuum optical depth is 2/3. Using the critical point, r c l , to test the solution 
for the second dissipation law given by v 2 produces a smaller mass-loss rate. At that choice 
of critical point, the velocity v c is larger than that predicted for v\ according to Equation 
[TBI Similarly the derivative of the velocity at the critical point is larger than in the case of 
the first dissipation law. At the star's surface, r = R*, the dissipation law does not play a 
role, implying the momentum equation and the quantity v(R*)v'(R*) is approximately the 
same for both cases. Hence there are three possibilities for the velocity at the surface: the 
velocity for the second dissipation law is less than, equal, or greater than the velocity for 
the first dissipation law. If the velocity given by v 2 is less, then the velocity derivative at 
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the surface is greater, which causes the velocity as a function of distance from the stellar 
surface to increase too rapidly and not satisfy the Singularity and Regularity conditions at 
r c . If the velocity at the surface for v 2 is greater than or equal to that in the first case then 
the density structure given by p oc M/v 2 is smaller at all distances. The continuum optical 
depth is less than 2/3 for the second dissipation law and thus the critical point r c that would 
satisfy the second dissipation is at a larger distance from the star than the critical point in 
the first case. The larger value of r c for the case of v 2 increases the mass-loss rate and acts 
to cancel the effect of greater dissipation. 

The conclusion of these two sections is the dissipation power law will have a significant 
effect on density and velocity structure but only a minor effect on the mass-loss rate. The 
effect on mass loss due to dissipation given by the power law is cancelled by forcing a different 
value of the critical point to satisfy the requirement that the continuum optical depth of the 
wind is 2/3. So far, the dissipation exponent v has been treated as a variable where specific 
cases have been explored, v = 0, 2, 3. The results show that u, while not greatly affecting 
the mass-loss rate, causes complications regarding the Parker point, r p , except for the case 
of v = 2. Therefore, to avoid those pitfalls, the mass loss of pulsating winds will be treated 
with the dissipation law given by v = 2. 



2.5. Comparison of Pulsating and Non— Pulsating Winds 

It is useful to compare the instantaneous mass-loss rate at some phase of pulsation to 
the mass-loss rate predicted by only radiative driving in a simple case. In this case, one 
can consider the region of the wind where v » a, ignoring terms of the order a 2 /v 2 and 
assuming the term 2a 2 r « GM»(1 — T e ) — (Rl (again v is assumed to equal 2). These 
assumptions still allow the solution to prob e most of the subcritical r egion r < r c as well 



as the supercritical region r > r c to large r (ILamers &: Cassinellil Il999h . but does limit the 



solution from describing the case where the GM*(1 — Y e ) — (R 2 « 0. Thus it is possible to 
rewrite Equation [TBI as 

W c « (r^A [gm,(i - r e ) - (Rl] ■ (46) 

If it is assumed that the right hand term, over the range specified, does not vary then r c v c v' c 
is approximately constant meaning rvv' is constant over that range. Therefore the velocity 
structure of the wind ejected at some instant can be approximated as 



2 2 
V =V C 



(_^_) [G M.(i-r e )-c«](i-i). (47) 
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Furthermore the term v 2 can be replaced with the approximate form of Equation [T5], such 
that the velocity SIS db function of r is 



I) 



ft 



1 — ft 



Rearranging and solving for r c at the surface r = i?* and taking v(R*) ~ a 

-l 



rv t 



(4$ 



1 + 



1 — ft a 



ft f. 



esc,efF 



(49) 



where vl sceS = GM»(1 — T e )/R* — CR*- This result is simila r to the solution of the critical 



radius for radiative driving from lLamers Cassinelfil (119991 ) . except the escape velocity is 
different. 

When the function ( is greater than zero then the effective escape velocity for the 
pulsating case is smaller than the escape velocity when there is no pulsation implying r c in 
the case of pulsating stars is smaller than the critical point for no pulsation. Therefore in 
the quasi-static pulsating limit one would expect the wind to have a larger acceleration than 
in the case of only radiative driving. The instantaneous mass-loss rate in the approximate 
case, using the assumptions given above with Equation [T5J is 



M 



a e L*k Z 
4vrc ~ZZ 



(1 



ft. 



l/a 



An 



a 



1 



a 



[GM*(l-r e )-{Rl] 



l-l/a 



(50) 



Mp U l s 


~<Je(ty 


-l+l/a 




l/a / \ 






l-l/o 


M rad 












[ GM,(i-r e )_ 





/ \VeVth, 

Relative to the approximate mass-loss rate due to radiative driving only, the ratio of mass- 
loss rates is 



(51) 



where the quantities denoted with a bar are time averaged over one period of pulsation. It 
is interesting to note the instantaneous mass-loss rate scales directly with the variations of 
luminosity to the power of 1/ft, and the electron scattering opacity to the power —l + l/a 
and inversely with the square root of the temperature through the thermal velocity. More 
important is the last term in square brackets containing the ratio of the acceleration due to 
shocks and pulsation and the effective gravitational acceleration. The exponent 1 — 1/ft is less 
than zero as a is of order 1/2, meaning as ( becomes larger so will the instantaneous mass- 
loss rate. Furthermore it is clear when ( is of order the effective gravity the mass-loss rate 
goes to infinity; however, this regime violates the assumption 2a 2 r << GM t (l — T e ) — (R 2 . 
Still this does illustrate the non-linear behavior of the mass-loss rate due to the added effects 
of pulsation and shocks. The goal now is to quantify how large the Cepheid mass-loss rates 
can be. 
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3. Testing the Wind Model with Cepheid Data 



The Cepheid instability strip spans a large range of mass, radius and luminosity, so to 
understand the role of mass loss, it is necessary to have a significant number of Cepheid 
models and observations. That information can be used to calculate mass-loss rates across 
the instability strip and probe how mass loss evolves with time. The circumstellar shells 
that have been observed can be modeled as material from the wind that condenses to dust 



when the temperature reaches approximately 1500 K, as mentioned by lKervella et all (120061 ). 
The dust would absorb light and re-emit isotropically, some of the radiation is toward the 
observer and thereby causing a flux excess at near infrared and longer wavelengths. 

To calculate mass-loss rates using the modified version of the CAK method, it is neces- 
sary to know the following parameters: the radius, amplitude of radius variation, luminosity, 
amplitude of luminosity variation, mass, and the period, along with the metallicity. Fur- 
thermore the values of a and k are needed for radiative driving. These parameters are based 
on the effective temperature and gravity of the star . For th i s wor k the parameters chosen 
are: a = 0.465 and k = 0.064 from the analysis of lAbbottl (Il982l ). These values are used 
for models with temperatures ranging from 4600 K to 6000 K. There is a lack of parame- 
ters calculated at effective temperatures less then 6000 K, but a and k are not likely to vary 
much since the main elements that drive the mass loss radiatively are ir on, neon and calcium 
as we ll as contributions from hydrogen and helium remain the same (ILamers fc Cassinelli 
19991 ). Molecules ma y play a small role as molecular opacities are important in AGB star 
winds, T e fj « 3000 K ([Helling et al.ll2000l ). but they are ignored in this analysis. 



The amplitude of the ra dius variation can be d eterm ined by integrating the radial ve- 
locity profile over one period. iMoskalik fc Gorynyal (120051 ) compiled a list of viable interfer- 
ometric targets and calculated the amplitudes of radius variation. The Cepheids on this list 
that pulsate in the fundamental mode are used here. Since all of the Cepheids in the sample 
are galactic stars the me t allicit y is assumed to be solar, Z = 0.02, consistent with the shock 
model from iFokin et all (119961 ) . Recently, lAsplund et all (120051 ) have suggested that the 
solar metallicity might be Z & = 0.012, This downward revision is being questioned because 
it is incompatible with helioseismology. It is unclear how the lower metallicity would affect 
the shock structure in the atmosphere of a Cepheid, although if one increases the hydrogen 
and helium abundances then the shocks generated in the hydrogen and helium ionization 
fronts might have more energy. An obvious consequence of a lower metallicity is seen in 
Equation [501 where the mass-loss rate explicitly depends on the metallicity, and would lead 
to the mass-loss rate at lower Z being 1/3 the mass-loss rate at Z = 0.02. However, it is 
shown other parameters affect the mass loss more significantly. The radius, luminosity and 
luminosity variation are found in the David Dunlap Observatory Catalogue of Galactic Clas- 
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Name 


Period(d) 


My 


AV 


Radius (R©) 


AR/R 


Mass (M©) 


Teff(K) 


R TrA 


3.388 


-2.74 


0.561 


25.3 


0.074 


2.4 


6470 


RT Aur 


3.724 


-2.86 


0.803 


35.1 


0.103 


4.7 


5640 


BF Oph 


4.064 


-2.97 


0.636 


34.5 


0.101 


3.9 


5840 


V Vel 


4.375 


-3.06 


0.689 


32.8 


0.116 


3.1 


6110 


T Vul 


4.436 


-3.08 


0.643 


38.2 


0.100 


4.5 


5690 


V482 Sco 


4.529 


-3.11 


0.652 


44.4 


0.103 


6.4 


5320 


S Cru 


4.688 


-3.15 


0.690 


42.1 


0.108 


5.3 


5510 


AP Sgr 


5.058 


-3.25 


0.832 


44.0 


0.110 


5.3 


5510 


V350 Sgr 


5.152 


-3.27 


0.705 


47.8 


0.110 


6.3 


5320 


6 Cep 


5.370 


-3.32 


0.838 


41.6 


0.116 


4.2 


5800 


V Cen 


5.495 


-3.35 


0.804 


45.3 


0.119 


5.0 


5820 


Y Sgr 


5.768 


-3.41 


0.725 


50.0 


0.120 


6.0 


5370 


RV Sco 


6.067 


-3.47 


0.824 


47.7 


0.102 


4.9 


5570 


S TrA 


6.324 


-3.53 


0.768 


39.2 


0.106 


2.8 


6230 


AW Per 


6.456 


-3.55 


0.812 


47.3 


0.118 


4.4 


5700 


BB Sgr 


6.637 


-3.59 


0.597 


40.3 


0.097 


2.8 


6230 


AT Pup 


6.668 


-3.59 


0.904 


45.4 


0.120 


3.8 


5870 


V Car 


6.699 


-3.60 


0.601 


40.6 


0.094 


2.8 


6220 


U Sgr 


6.745 


-3.61 


0.717 


51.4 


0.115 


5.1 


5540 


V496 Aql 


6.808 


-3.62 


0.349 


45.5 


0.067 


3.7 


5900 


X Sgr 


7.014 


-3.66 


0.590 


49.8 


0.093 


4.4 


5720 


U Aql 


7.031 


-3.66 


0.757 


54.7 


0.116 


5.6 


5440 


77 Aql 


7.178 


-3.69 


0.799 


54.9 


0.121 


5.5 


5510 


W Sgr 


7.603 


-3.76 


0.805 


63.3 


0.121 


7.2 


5200 


RX Cam 


7.907 


-3.81 


0.729 


76.0 


0.120 


10.8 


4770 


W Gem 


7.907 


-3.81 


0.822 


50.7 


0.126 


3.9 


5840 


U Vul 


7.998 


-3.82 


0.718 


56.5 


0.098 


5.0 


5550 


GH Lup 


9.268 


-4.01 


0.192 


58.2 


0.040 


4.3 


5710 


S Mus 


9.660 


-4.06 


0.500 


71.3 


0.109 


6.8 


5220 


S Nor 


9.750 


-4.07 


0.640 


66.4 


0.119 


5.6 


5420 


(3 Dor 


9.840 


-4.08 


0.630 


64.4 


0.118 


5.1 


5550 


C Gem 


10.139 


-4.12 


0.480 


64.9 


0.099 


5.0 


5590 


XX Cen 


10.965 


-4.22 


0.924 


57.8 


0.140 


3.3 


6020 


RX Aur 


11.614 


-4.29 


0.664 


63.4 


0.122 


3.8 


5840 


TT Aql 


13.740 


-4.51 


1.082 


95.8 


0.196 


8.5 


5020 


X Cve 


16.368 


-4.73 


0.986 


118.1 


0.216 


11.1 


4760 


Y Oph 


17.139 


-4.78 


0.483 


93.5 


0.080 


5.8 


5410 


VY Car 


19.011 


-4.91 


1.065 


108.9 


0.220 


7.3 


5140 


RZ Vel 


20.417 


-5.00 


1.181 


111.7 


0.243 


7.0 


5210 


T Mon 


27.040 


-5.36 


1.028 


130.6 


0.230 


6.8 


5230 


1 Car 


35.563 


-5.70 


0.725 


183.7 


0.191 


10.7 


4780 


U Car 


38.726 


- 5.81 


1.165 


162.2 


0.238 


6.9 


5200 


RS Pup 


41.400 


-5.90 


1.105 


197.7 


0.246 


10.3 


4820 


SV Vul 


44.978 


-6.00 


1.054 


235.5 


0.246 


14.2 


4520 



Table 2: Data for modeling the mass-loss rate behavior of Cepheids. See text for description 
and references. 



- 21 - 



Name 


M puls {M Q /yr) 


M rad 


{M Q /yr) 


Name 


M puls (M e /yr) 


M rad (M /|/r) 


R TrA 


4.9 x 10- 


-10 


1 


2 


x 


10' 


-ii 


rj Aql 


1.8 x 10- 


10 


3.6 x 


io- 


-11 


RT Aur 


2.3 x 10 


-8 


7 


5 


x 


10" 


-12 


W Sgr 


1.3 x 10- 


10 


3.0 x 


io- 


-11 


BF Oph 


5.4 x 10" 


-10 


1 


1 


X 


10" 


-11 


RX Cam 


8.4 x 10" 


11 


2.0 x 


io- 


-11 


V Vel 


1.3 x 10 


-8 


1 


7 


X 


10" 


-11 


W Gem 


3.0 x 10- 


10 


6.2 x 


io- 


-11 


T Vul 


1.3 x 10" 


-10 


1 


2 


X 


10" 


-11 


U Vul 


1.3 x 10- 


10 


4.8 x 


io- 


-11 


V482 Sco 


1.0 x 10" 


-10 


9 





X 


10" 


-12 


GH Lup 


1.4 x 10- 


10 


8.3 x 


io- 


-11 


S Cru 


2.1 x 10- 


-10 


1 


2 


X 


10" 


-11 


S Mus 


1.4 x 10- 


10 


5.6 x 


io- 


-11 


AP Sgr 


1.5 x 10- 


-10 


1 


4 


X 


10" 


-11 


S Nor 


2.0 x 10" 


10 


7.0 x 


io- 


-11 


V350 Sgr 


1.2 x 10- 


-10 


1 


2 


X 


10" 


-11 


/3 Dor 


2.4 x 10- 


10 


8.3 x 


io- 


-11 


5 Cep 


3.6 x 10- 


-10 


2 


3 


X 


10" 


-11 


C Gem 


2.2 x 10" 


10 


9.5 x 


io- 


-11 


V Cen 


1.7 x 10 


-9 


2 


7 


X 


10" 


-11 


XX Cen 


7.3 x 10" 


10 


1.7 x 


io- 


-10 


Y Sgr 


1.9 x 10 


10 


1 


7 


X 


10" 


-11 


RX Aur 


4.6 x 10- 


10 


1.7 x 


io- 


-10 


RV Sco 


9.3 x 10- 


-11 


2 


4 


X 


10" 


-11 


TT Aql 


8.0 x 10' 


-9 


1.1 X 


io- 


-10 


S TrA 


2.3 x 10- 


-10 


5 





X 


10" 


-11 


XCyg 


7.7 x 10- 


-8 


1.3 x 


io- 


-10 


AW Per 


1.8 x 10- 


-10 


3 


2 


X 


10" 


-11 


Y Oph 


5.1 x 10- 


10 


3.0 x 


io- 


-10 


BB Sgr 


2.1 x 10- 


-10 


5 


7 


X 


10" 


-11 


VY Car 


4.7 x 10" 


-9 


2.8 x 


io- 


-10 


AT Pup 


2.4 x 10- 


-10 


1 


1 


X 


10" 


-11 


RZ Vel 


3.7 x 10- 


-8 


3.8 x 


io- 


-10 


V Car 


2.0 x 10- 


-10 


5 


8 


X 


10" 


-11 


T Mon 


5.5 x 10- 


-9 


8.1 x 


10 


10 


V496 Aql 


1.0 x 10- 


-10 


4 


5 


X 


10" 


-11 


I Car 


2.4 x 10- 


-9 


9.9 x 10- 


-10 


X Sgr 


1.2 x 10- 


-10 


4 


1 


X 


10" 


-11 


U Car 


1.0 x 10- 


-8 


2.0 x 10 


-9 


U Sgr 


1.4 x 10- 


-10 


3 


1 


X 


10" 


-11 


RS Pup 


6.5 x 10- 


-9 


1.6 x 10 


-9 


U Aql 


1.3 x 10- 


-10 


3 


1 


X 


10" 


-11 


SV Vul 


4.8 x 10- 


-9 


1.3 x 10 


-9 



Table 3: Predicted mass-loss rates for Cepheids using shock and pulsation dynamics along 
radiative driving and assuming only radiative driving. 



sical Cepheids (IFernie et al.lll995l ). The effective temperature iscalculated from the mean 
lumi nosity and rad ius while the mass is ca lculated using the Period-Mass-Radius relation 
fromlGiererJ Jl989h and lFricke et~aD Jl972h : 



M* = [4OP(days)]- L49 ( J R,/ J R ) 2 - 53 M . 



\2.53 



(52) 



The properties of the observed Cepheids are listed in Table [2] and the computed mass-loss 
rates are listed in Table [3J 

The first test is to compare the sensitivity of the pulsation mass loss to the value of 
the critical point, r c . If the mass loss varies significantly as a function of r c for a given set 
of parameters then the predicted mass-loss rate will strongly depend on the criterion that 
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Fig. 1. — (Left Panel) The predicted continuum optical depth of the wind for S Cep and / 
Car at minimum radius. The horizontal lines highlight the regime where the r is of order 
2/3. (Right Panel) The predicted mass-loss rate for 5 Cep and I Car at minimum radius. 
The rates do not vary significantly with r c . 



the continuum optical depth of the wind be 2/3. The continuum optical depth is shown in 
Figured] (Left) as a function of the critical point, r c , for two Cepheids in the sample, I Car 
and <5 Cep at minimum radius. The continuum optical depth varies between 0.1 and 1 over 
a small range of critical points, for 5 Cep the range is approximately 0.5-R* and for I Car it 
is about 1.5-R*. Over the two ranges of critical points, the mass-loss rates vary only about a 
factor of 2 and 3 as shown in Figured] (Right). Therefore, the predicted mass-loss rate of a 
Cepheid is not sensitive to the criterion requiring the continuum optical depth of the wind 
to be 2/3. 

As a reference, mass-loss rates for this set of observational data are predicted assuming 
a static stellar surface and the only driving mechanism is radiative with pulsation terms 
ignored. In Figure [2J the mass-loss rates are plotted as a function of pulsation period. The 
linear relation in the Log-Log plot is striking, if one considers the equation for mass loss in 
the CAK method in the approximate case given in Section 2 and ignoring pulsation, then 

logM « ( 1 ] logM* + -logL - -logT eff + Const. (53) 

\ a J a 2 

The mass can be written in terms of the period and radius via Equation [52] and the radius 
is a function of the effective temperature and luminosity, R oc L l l 2 T~ 2 . Furthermore the 
effective temp e ratur e can be represented by the color (B — V)q using the transform from 



Fry fc Carney! (]1999l ). The radiative driven mass loss is thus 



log M rad (M & /yr) = -14.047 + 0.695 log L(L Q ) - 0.574(5 - V) + 1.71 log P{d). (54) 
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Fig. 2. — The predicted mass-loss rates for Cepheids in the quasi-static approximation usin^ 
radiative driving as the only mechanism for generating the wind. 



The lu minosity can be eliminated using a Period-Luminosity-Color relation from ldi Benedetto 
(I1995I ) to obtain an MPC relation for the radiative driven mass-loss rates for specific 
Cepheids, 

logM rad = -12.11 + 2.69 log P - 1.16(5 - V) . (55) 

This relation is for specific Cepheids but to understand the apparent relation of radiative 
mass loss as a function of period in Figure we wish to eliminate the color dependence 
to derive a statistical representation of the mass-loss rate. Returning to Equation [5U the 
luminosity an d the color can be elim inated by using Period-Luminosity and Period-Color 
relations from iTammann et al.l (120031 ) yielding 



logM rad = -12.45 + 2.27 log P. 



(56) 



This relation is a theoretical derivation of the radiative driven mass-loss rates for the sta- 
tistical sample of galactic Cepheids based on best-fit Period-Luminosity and Period-Color 
relations. 

Performing a least squares fit on the predicted radiative driven mass-loss rates as a 
function of the period as shown in Figure [H one finds 



log M rad (M 0/ V) = 2.11ogP(d) - 12.2. 



(57) 



Fig. 3. — (Left Panel) The mass-loss rates for Cepheids using the combination of radiative 
driving and pulsation but ignoring shocks to generate the wind plotted with the mass-loss 
rates from only radiative driving for comparison. (Right Panel) The ratio of the mass-loss 
rates computed using pulsation and that using only radiative driving. Accelerating the wind 
via pulsation does not effect the mass loss significantly by < 50%. The mass loss is more 
enhanced for Cepheids with period > 10 days is due to the combination of larger amplitudes 
of radius variation on average and lower gravity. 



The two relations differ but this is to be expected. The derived result ignores the contri- 
butions due to continuum radiative driving and assumes that GM* » a 2 r. The best fit 
relation will have statistical uncertainties as the majority of the Cepheids examined have 
periods P < lOd. Therefore the best-fit and derived relations for the radiative driven mass- 
loss rate as a function of period agree and provide a rough lower limit for the mass-loss 
rates of Cepheid s. Furthermor e this best-fit relation is equivalent to a Reimer's relation 



M = 1]L*R*/M* ( lReimerslll977l ). where all quantities are in solar units and time is in years. 
The standard value is r\ = 10~ 13 ; for Cepheids this would predict mass-loss rates much 
larger than found here using the CAK method. Fitting this relation to the radiative driving 
mass-loss rates determined here, one finds a value of 77 ~ 4.4 x 10~ 15 . 

The second step is to compute the mass-loss rates of Cepheids assuming the contri- 
bution from shocks is zero and the wind is accelerated by momentum from pulsation and 
radiation. The result is shown in Figure [3] (Left) where the mass-loss rates due to pulsation 
but not shocks are plotted alongside the mass-loss rates from radiative driving only. The 
comparison shows the pulsation does not greatly enhance the mass loss and the result is fur- 
ther highlighted in Figure [3] (Right) where the ratio of the two mass-loss rates is also shown. 
Pulsation does amplify the mass-loss rates but it is a 50% effect at most, which will not 
produce significant mass loss. The result is not surprising if one considers Equation [51] and 
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Fig. 4. — (Left Panel) The mass-loss rates for Cepheids using the combination of radiative 
driving, pulsation and shocks to generate the wind. The dotted line is the least square fit to 
the radiative driven mass-loss rates shown in Figure [2J (Right Panel) A comparison of the 
mass-loss rates computed using pulsation plus shock effects and that using only radiative 
driving. The wind is strongly enhanced for some Cepheids with the largest enhancement 
being a factor of approximately 750 and the smallest about 1.7 times. 



replaces ( with just u 2 AR. Cepheid pulsation can be approximated as a linear perturbation 
implying the acceleration due to pulsation is much smaller than the gravity of the Cepheid. 
In that case, the ratio of pulsation mass-loss rates to the radiatively driven mass-loss rates 
is close to unity. Therefore it is necessary to consider the effect of shocks on the wind. 

When the terms describing the shocks and pulsation are included in the calculation, the 
mass-loss rates change dramatically. This can be seen in Figure H] (Left) where the pulsation 
mass loss is plotted as a function of period. The mass-loss rates do not follow a simple 
relation such as a M — P relation for radiatively driven mass loss or a Reimer's relation. 
Some Cepheids have very large mass-loss rates compared to other Cepheids, which can be 
understood from the approximate version of the mass-loss rate, Equation [521 The enhanced 
mass-loss rate, shown in Figure H] (Right), depends on (1 — C/sO -1 if one assumes a ~ 1/2, 
where g is the gravity of the star; the term is maximized when either ( is large or g is 
small. The gravity is smallest when the Cepheid is at maximum expansion, = 0.5, but the 
acceleration due to pulsation is at a minimum and the shock amplitude is zero. Therefore ( 
at this phase, is acting to decrease the mass-loss rate. The pulsation plus shock function, (, 
is largest at minimum radius where the shock amplitude is large and the acceleration due to 
pulsation is maximum, meaning both effects contribute. The sum of the two accelerations, 
which act at similar magnitudes, and the fact that the pulsation depends on the amplitude 
AR implies there is no simple formula for the mass-loss rate. 
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Fig. 5. — The observed Cepheids plotted on a Cepheid instability strip of the HR diagram 
where the dashed lines represent approximate boundaries. (Left Panel) The size of the circles 
represent the calculated mass-loss rates for the Cepheids and (Right Panel) the size of the 
squares represent the ratio of the pulsational mass loss and the radiative driven mass loss. 
All of the brighter Cepheids, log L/L Q > 3.7, have large mass-loss rates, while this is true 
for only a few of the less bright Cepheids but only a fraction of the luminous Cepheids have 
their mass loss significantly increased due to pulsation and shocks. 



Only some Cepheids have a dramatic enhancement of mass loss, implying there are 
restricted regions in the instability strip where Cepheids are more susceptible to lose mass. 
Figure [5] (Left) shows the location of the Cepheids modeled here with the size of the symbols 
representing the amount of mass loss. The plot shows that larger Cepheids existing in the 
upper part of instability strip exhibit higher mass-loss rates, as would be expected from the 
calculations of radiative driving. It is also interesting that the few short period Cepheids with 
high mass-loss rates are scattered along the effective temperature axis but have consistently 
lower luminosity, which could be a result of a lower mass and short period contributing 
to lower gravity and higher pulsation plus shock acceleration respectively. Figure [5] (Right) 
shows the Cepheids on the HR diagram but with the point sizes now representing the amount 
of mass loss enhancement. The results are striking; there appears to be two bands where 
the mass loss is enhanced. This may imply the mass loss may be related to which crossing 
of the instability strip the Cepheid is making. 

There are many sources of uncertainty in this analysis. The most important source 
is the mass for each Cepheid, based on a Period-Mass-Radius relation, which has scatter. 
Deviations of the mass in the relation to the observed values are i mportant. Conside r the 



Cepheid S Mus, for which a mass of 6.8M is used here; however lEvans et al.l (120061 ) find 



S Mus has a mass of 6.0 ± O.4M . It is clear from Equation [51] that a lower mass will 
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Fig. 6. — The predicted mass-loss rate for S Mus found by holding all necessary parameters 
describing the Cepheid constant except for the mass, which is allowed to vary. 



increase the rate of mass loss if all other parameters remain the same. One can test the 
dependence of the mass-loss rate on the mass of the Cepheid by calculating a set of models 
with varying mass and holding all other parameters the same. The effect is shown in Figure 
[6] for S Mus, ranging the mass from 5.6M to 6.8M . For the larger masses the mass-loss 
rate is not sensitive to the mass, at 6.8M , M = 1.4 x 10~ 10 M Q /yr and increases by only a 
factor of five at 6.OM . At lower masses, however, the mass-loss rate exhibits a non-linear 
dependence on the mass increasing by two orders of magnitude from M = 6.0 to M — 5.6M . 
Therefore when the mass-loss rates are large, it is important to have strong constraints on 
the mass. 

The uncertainty in the calculation may be quantified by considering the analytic deriva- 
tion of the error of the pulsation-driven mass-loss rate based on Equation [THJ simplified 
by assuming the quantities a 2 r c and T e are insignificant. Writing the error as AF(xi) = 
y^JdF '/dxi) 2 al t where a x , is the error of each quantity, 

M 

In this equation, it is also assumed that the uncertainty of the period and amplitude of 



aL, 



2T, 



G 2 a 2 M + Ria 2 + 4( 2 Kia 



I 



R 



(gm* - cm 2 



(58) 
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luminosity variation is negligible. The uncertainty of the function ( is given by 



cr f = <J [uj 2 a AR \ 2 + 



du$ 
Pd4> 



L 



1/8 



P 

Ps 



-7/24 



(-) 



5/12' 



2 r 



-ZIl\ + ( 5 17 M 
8LJ V 12 M* 



1/2 



(59) 

Upon considering the third term of Equation EE, the uncertainty of the mass-loss rate is 
inversely dependent on the balance of forces GM t — (Rl. It was shown in Section 2 that 
the pulsation driven mass-loss rate is large when the balance of forces is small and likewise 
a A -j/M is also inversely proportional to the balance of forces. This implies that when the 
mass-loss rate is large due to pulsation and shocks then the error is even larger. In this 
case, the uncertainty may be larger by more than an order of magnitude depending on 
the uncertainty of the mass, radius, luminosity, and amplitude of radius variation. The 
uncertainty in the effective temperature is not a significant source of error. This implies that 
the calculation is very sensitive to the values of the parameters used and reiterates the need 
for more precise values of mass, as well as luminosity and radius. 

In this work, the uncertainty of the luminosity is about 10%, the uncertainty of the 
mass is about 25% to account for the error in the PMR relation, and the uncertainty of 
the am plitude of variation o f radius is about 5% based on the uncertainty of the projection 
factor (iNardetto et al.l 120071 ). The uncertainty of the effective temperature is ±200f^ while 
the uncertainty of the radius is calculated from the luminosity and effective temperature 
errors. The fractional error of the pulsation driven mass-loss rate is shown in Figured (Left) 
as a function of period and the mass-loss rates are plotted as a function of period with error 
bars in the right panel. This shows that the error of the mass loss is sensitive to the balance 
of gravity, and the uncertainty of the mass and luminosity has little effect unless they are 
at least an order of magnitude less. This reiterates the importance of the balance of forces 
in determining the mass-loss rates. While the uncertainties of the predicted values of the 
mass-loss rates are large in some cases, the observational determinations of mass-loss rates 
(as shown in the next section) also have large uncertainties. 



4. Comparison to Observed Mass— Loss Rates 

There have been only a small number of measurements of mass-loss rates for Cepheids. 
There are two main methods: observing the near and mid-infrared flux excess from dust 
and using that to determine mass loss, or from emission lines inferring large velocities and 
particle densities to measure mass loss. The Cepheids with estimated mass-loss rates that 
are coincident with those used in this work are listed in Table HJ By comparing the measure- 
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Fig. 7. — (Left) The fractional errors of the pulsation driven mass-loss rates for the sample 
of galactic Cepheids as a function of period. (Right) The mass-loss rates of the Cepheids 
as a function of period with errors calculated using Equation [58] added. The errors appear 
related to the enhancement of the mass loss due to pulsation and shocks and the plot shows 
that uncertainty of the mass-loss rate has a similar non-linear behavior as a function of 
mass as the mass-loss rate which is shown in Figure El 



Name 


M(M Q /yr) 


Reference 


RS Pup 


3.5 x KT 6 


Deasv (1988) 


/ Car 


2 x KT 8 


Bohm-Vitense & Love 1994) 


R TrA 


3 x KT 9 


McAlarv & Welch (1986) 


S Mus 


< 10- 9 


Rodrieues & Bohm-Vitense (1992) 


5 Cep 


< 5.5 x 10- 9 


Welch & Duric 1988) 


f] Aql 


< 6.1 x 10-9 


Welch & Duric 1988) 


T Mon 


< 4.2 x KT 8 


Welch & Duric 1988) 



Table 4: The inferred mass-loss rates of various Cepheids from the literature. 



ments with the predicted mass-loss rates including pulsation and shocks, it is apparent the 
predictions are lower than the measured rates by about a factor of 10, except for RS Pup 
which is a several orders of magnitude different. 



For the case of I Car, iBohm-Vitense &: Love! (11994 ) used ultraviolet spectra to detect 
emission lines of C II, C IV, Mg II and O I. They argued carbon emission lines provide evi- 
dence for mass loss because they require velocities of order 100 km/s to form, which is the 
same order of magnitude as the escape velocity. The Mg II lines are seen to have two emission 
components surrounding a broad central absorption profile that is relatively constant with 
pulsation phase. The width of the central absorption profile of magne sium lines is believed to 
be a result of circumstellar material or interstellar material or both. Bohm-Vitense fc Love 
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( 119941 ) use the velocities inferred from the lines to calculate an optical depth for the ab- 
sorption shell and thus a column density. The column density depends on the choice of the 
turbulent velocity, and the uncertainty in the turbulent velocity causes an exponentia l change 
in the estimate of the column density; for this reason Bohm-Vitense fc Lovel (1994) argued 
the m ass-loss rate was uncertain by two orders of magnitude. Because iBohm-Vitense fc Love 
(11994] ) used a value for the turbulent velocity near the lower limit, their mass-loss rate is an 
upper limit. The mass loss derived by these authors depends also on the dis tance to I Car, 
which was taken to be d = 400 pc. However, recent parallax measurements ([Benedict et al. 



20071 ) suggest I Car is about 500 pc away, thus lowering this upper limit as well. Given these 
qualifications, the value of the mass-loss rate p redicted here, 2.4 x 10~ 8 M & /yr, is consistent 
with the result of IBohm-Vitense &: Lovel ( 19941) . 



IUE observations of S Mus (IRodrigues fc Bohm-Vitensd Il992l ) showed features in the 
absorption profile due to the Cepheid, its B5V companion and possibly the wind from the 
Cepheid. The authors modeled the line profile by assuming a velocity law for the wind, 
which they assumed to be spherically symmetric and unperturbed by either pulsation or 
the companion. Using a /3-law for the wind, t>oo(l — R*/r)^, gives a mass-loss rate of 
approximately 1.1 — 1.3 x 10~ 10 M Q /yr by fitting the observed absorption profile. Using 
different velocity laws, such as exponential and power-laws, had a dramatic effect on the 
mass-loss rate, with the maximum estimate being 2.5 x 10 5 times larger than the minimum 
estimate. The result from using the (3— law alone is uncertain by a factor of about 6 according 
to the authors, due to issues of fitting the continuum and the assumed abundances for the 
system. Therefore the result is an optimistic upper limit. The important conclusion is that 
our predicted mass loss is small like the observed value, tho ugh both may be wrong if the 
lower limits of the mass determined by lEvans et al.l (120061 ) is the true value of the mass 
as shown in Figure El The mass- loss rate will still be consistent with that calculated by 
Rodrigues fc Bohm-Vitensd (119921 ) . 



McAlary fc Welch! (119861 ) used IRAS observations to detect Cepheids and found modest 



infrared excess in several. By assuming the infrared excess is due to dust that formed in a 
stellar wind, the authors estimated mass-loss rates of the order 10~ 9 — 10~ 8 M & /yr and worked 
out the case of R TrA. The calculated mass-loss rate for R TrA is 4 x 10 -9 M Q / yr , about a 
factor of 9 greater than that found here. One of the differences in the m ass-loss calculation 
may b e due to differences in the choice of stellar radius and temperature, iMcAlary fc Welch 
( 119861 ) used a larger radius and smaller temperature. A larger radius will affect the calculated 
density of the shell while the smaller temperature will predict a larger flux excess. Therefore 
the two estimates for the mass-loss rate of R TrA are probably consistent. 



Analysis of IRAS observations was repeated by iDeasyl (119881 ) who found mass-loss rates 
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consistent with the previous work. Furthermore iDeasyl (119881 ) derived a mass-loss ra te from 
the ge ometric structure of the surrounding nebula for RS Pup. The observations of lHavlen 
( 119721 ) showed the reflection nebula was associated with the Cepheid and modeled the nebula 
as a set of four concentric dust shells, each having a mass of 0.05-0. 1M . By assuming a 
velocity for the stellar wind, it is possible to determine a timescale for the shells, but the 
timescale of the innermost shell is much longer than for the others. This implies the existence 
of an undetected shell based on the argument the shell timescales are related to the crossing 
time of a Cepheid on the instability strip. The mass-loss rate is approximated by the dust 
mass of the shell and the lifetime of the crossing of the instability strip. By assuming a shell 
mass based on the spacing of shells, the mass of the innermost shell may be overestimated, 
especially as it has not been detected and the shell is still being generated. Therefore the 
estimate is a measure of the upper limit of the mass-loss rate. 



Welch fc Durid (119881 ) observed a sample of Cepheids at radio wavelengths to measure 
mass-loss rates. They modeled an ionized wind at temperatures of order 10 4 -10 5 K. The 
mass-loss rate is not sensitive to this assumption, but the temperature should not exceed 
10 4 K, the temperature where shocks would be generated in a Cepheid. The authors also 
assume a wind velocity of order 100 km/s, which is not necessarily true. The value is 
chosen because the observed components of the wind must exist at a radius greater than 
4R*, meaning a lower wind velocity can be used if the components are older and have moved 
a further distance from the star. There is also systematic uncertainty about the choice of 
distances to these Cephe i ds; parallax measurements indicate a difference of order 20-30% 



(IFeast fc Catchpold 119971 ; Benedict et al 



2002 



although the model devised by Welch &: Durid ( 1988 



2007 ). The mass-loss rates are consistent, 
) is different than the one used here. 



The main conclusion to be drawn from these observed mass-loss rates is they are upper 
limits, with large uncertainties. These upper limits are orders of magnitude larger than the 
predicted mass-loss rates derived by assuming only radiative-driving, implying there must 
be additional driving forces, lending credence to the pulsation + shock mechanism described 
in this work. 



5. Model Infrared Excess 

The works of iMc Alary fc Welch! J 1986^ andlDeasv^ (119881) bo th discuss the i nfrared flux 
excess of Cepheids in the IRAS bands. iKervella et all tood) andlMerand et all (bood . l2007h 
discovered excess K-band flux using interferometry, and lEvans et al.l (120071 ) discuss using 
the Spitzer Space Telescope to search for IR excess. Predictions of infrared excess provide a 
useful test of the mechanism for mass loss in Cepheids. 
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Infrared excess can be produced in two distinct ways: by hot ionized winds or by cool 
dusty winds. Since the wind model in this work assumes it is driven by mechanical energy, 
the wind is not hot so a dusty wind is presumed. The wind is in radiative equilibrium with 
the Cepheid, so the temperature is dependent on the stellar temperature and is inversely 
proportional to the square of the distance from the star. Dust forms in the wind when the 
temp erature is near 1500 K, the condensation temperature for dust. Following the discussion 
from lLamers fc Cassinelfil (119991 ) for the luminosity of an optically thin dusty wind, one can 
represent the temperature structure as 

/ 2r x ^ 5 



T{r) 



Teff u 



(60) 



By rearranging, one can determine the condensation radius: for S Cep this is about 14.7/2*, 
while for the coolest Cepheid in the sample, SV Vul, the condensation radius is about 7.9/2*. 

To calculate the dust luminosity, it is assumed that the dust is forming far enough from 
the star that pulsation and shocks will not affect the structure. From the above calculation, 
however, this is not an ideal assumption, at least in the inne r edge of the dust sh e ll. T he 
total luminosity of the dust at a given frequency is derived by lLamers fc Cassinellil (Il999l ). 



3<o 2 >lM„ 



Qt jT" B u (T d ) jl - \ (l - y/1 - {Rjry) I 



dr. 



47r < cr > p d v d ./,. 12' 

where the quantities with a subscript d refer to dust. The terms in this expression are 
evaluated in the following way. The dust is assumed to be graphite with pd = 2.2g/cm 3 . 
The mass-loss rate of dust is one-hundredth of the total mass-loss rate, based on the galactic 
gas-to-dust mass ratio. Inside the condensation radius the value of is zero. The velocity, 
Vd, is the velocity of the wind at the distances from the star being considered. The value of Vd 
that is used is the mean terminal velocity of the wind averaged over one period of pulsation. 
To simplify the integration, r min = /?* and r max = oo; if the velocity is 100 km/s then the 
wind would travel a distance of order 100 pc in one million years. The dust will contribute 
to the infrared excess only at a much smaller distance, so th e assumption is reaso nable. The 
dust is assumed to have a grain size distribution as given by lMathis et al.l (119771 ). where the 
number density is n(a)da oc Ka~ 3,5 da. If the grain size ranges from 0.005 pm to 0.25 pm, 
then the term < a 2 > / < a 3 > is ~ 40 prrT 1 . The absorption efficiency, Q„ , is based on 
the argu ment that most of the dust absorption will be at optical wavelengths and will be of 
order 2 ( j Jones fc Merrill! 1 19761 ). Using these values, Equation [61] can be used to predict the 
luminosity of the dust shell in the infrared wavelengths of the VLTI, CHARA, IRAS and 
Spitzer. 

The infrared excess found from interferometric observations of Cepheids is summarized 
by lMerand et al.l (120071 ). where the authors list the fraction of the total flux contributed by 
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Fig. 8. — The predicted fraction of the total luminosity of a Cepheid contributed by the 
presence of dust at 2.2 fim. The circle represents the example of the smallest predicted 
circumstellar emission (CSE) flux, RX Cam, and the triangles represent RT Aur, V Vel, X 
Cyg and RZ Vel which have the largest predict CSE fluxes of the sample. 



circumstellar shells, defined as the circumstellar emission (CSE) flux. The results for the 
fundamental pulsating Cepheids are shown in Table [5j The non-pulsating supergiant a Per 
was observed as well and was found to have almost no contribution to the 2.2 fim flux due 
to circumstellar material. The infrared excess due to the predicted mass loss for the sample 
of Cepheids at 2.2 fim is shown in Figure El and the predicted CSE fluxes for the Cepheids 
observed with interferometry are given in Table [5] for comparison. These predictions are of 
the same order of magnitude as the observed, and differences may be due to the parameters 
in the dust model. The CSE flux is linearly dependent on the ratio < a 2 > / < a? >, 
the mass-loss rate of dust, terminal velocity, grain density, etc.; a small change of any of 
these parameters will change the prediction. The predicted CSE flux is most different for 
the cases of Y Oph and I Car, which could be related to the uncertainty of the inner radius 
of the dust shell or to the predicted mass-loss rates being underestimated for these two 
Cepheids; an increase of the mass-loss rate by a factor of 4 and 2 respectively would match 
the observations. In the enhanced CAK method this would require a decrease of chosen 
stellar mass of order of 1 to 1.5M Q or an increase of stellar radius of order 7-15i? , or some 
combination of the two, both of which are within the range of observational uncertainty. All 
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Name 


Observed CSE % 


Predicted CSE % 


5 Cep 


1.5 ±0.4 


1.3 


I Car 


4.2 ±0.2 


1.9 


Y Oph 


5.0 ±2.0 


0.95 



Table 5: The flux of circumstellar shells of Cepheids relative to the total flux observed using 
K-band interferometry compared to the predicted flux of circumstellar shells. 

things considered, the dust model is a consistent fit to the observations and the difference 
between the predicted and observed CSE is small. The predicted CSE fluxes of the sample 
Cepheids span a significant range, 0% to 50%, due to the combination of the fundamental 
parameters. This prediction can be tested. 

Th e IRA S observations provide fluxes in its four bands for a large sample of Cepheids. 



Deasyl (119881 ) found the ratio of F(25 fim)/F(12 fim) ranges from roughly 0.25 to 0.8. The 
lower end is the expected ratio from just blackbody radiation with no excess radiation (Figure 
1 of that paper). The ratio of the infrared excesses from the predicted mass-loss rates is 
given in Figure [9] (L eft). It is striking that the ratio of fluxes exhibit the same range of values 



as in [DeasyJ (119881 ). This agreement is encouraging as the luminosity of the dust at these 
two wavelengths is much less sensitive to the dust parameters, and the contributing dust 
is much farther from the inner boundary of the dust shell. The infrared excess predicted 
for just radiative driven mass loss is shown in Figured] (Right), where the maximum excess 
is significantly less than that predicted for pulsation plus shock driven mass l oss, a nd this 



prediction does not agree with the range of infrared excesses found by iDeasyl (119881 ) . This 
is further evidence that mass loss in Cepheids is driven by pulsation and not just radiation. 
Note there appears to be a minimum in the flux ratio in Figure [9] (Left) near logP ~ 0.9, 
which corresponds to log L/L Q ~ 3.5; this reflects the lack of mass-loss enhancement due to 
pulsation at that luminosity, which is consistent with the behavior as shown in Figure [5J 

Figures [8] and [9] (Left) show that the mass-loss rates and dust model can reproduce 
the existing interferometric and IRAS observations, but the Spitzer Space Telescope can 
now provide new results. Therefore it is useful to predict the infrared excess at wavelengths 
observed by Spitzer. The luminosity of the Cepheids with dust shells are plotted in Figure 
[TUl relative to the 3.6 /im luminosity, at wavelengths 8, 70, and 160 fim, respectively. The 
combination of the four wavelengths provide a test of mass-loss enhancement. At 8 fim, 
the excess is due to mass loss being enhanced by both pulsation and shocks and radiative 
driving at longer periods. This would provide a test for the period dependence of mass loss. 
Observations at the longer wavelengths provide a measure of smaller mass-loss rates because 
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Fig. 9. — (Left Panel) The ratio of the luminosity of the wind at 25 /im and 12 \im based on 
the predicted pulsation mass-loss rates. The ratio of the luminosity at these wavelengths is 
about 0.26 when there is no infrared excess, represented by the dashed line. The Cepheid, 
RX Cam, is represented by the circle at logP ~ 0.9 with the smallest infrared excess of the 
sample. The triangles represent the largest predicted flux excess for the Cepheids for RT 
Aur, V Vel, X Cyg and RZ Vel. (Right Panel) The ratio of the luminosity of the wind at 
25 /im and 12 /im for the predicted mass-loss rates based on radiative driving alone. The 
ratio of the luminosity at these wavelengths is about 0.26 when there is no infrared excess, 
represented by the dashed line. 



fluxes at these wavelengths are in the tail end of the stellar blackbody function as well as a 
measure of the size of the circumstellar shell. 

There are extreme values of the predicted flux excess; from Figures [9] (Left) and [10] it 
is clear one Cepheid, at period logP = 0.9, has minimal flux excess shown as a circle in the 
figures. This is the model for RX Cam, where the predicted mass-loss rate is 8.4 x 10 -11 , 
which is most likely due to the estimate of the mass M = 1O.8M . The infrared excess 
and mass-loss rate may easily be underestimated, but the result does show it is possible for 
Cepheids to have a low infrared excess similar to non-pulsating yellow supergiants. On the 
other hand, there are a number of Cepheids predicted to have very large IR excess, where the 
dust luminosity contributes most of the total luminosity at that wavelength. From Figure 
M (Left), there are four Cepheids with a CSE flux of order 50%: the short-period Cepheids 
RT Aur and V Vel and the 10-20 day period Cepheids X Cyg and RZ Vel. The two short 
period Cepheids have a large enhancement of mass loss, both by almost 10 3 times and hence 
the infrared excess is large. The same is true for the longer period Cepheids but to a lesser 
extent. It is unlikely the excesses are overestimated, if it is assumed the mass-loss rates 
are reasonable. Both V Vel and RT Aur have large effective temperatures for Cepheids, 
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Fig. 10. — The ratio of the predicted luminosity at 8.0 fim (Left Panel), 70 fim (Middle 
Panel), and 160 jim (Right Panel) relative to the predicted 3.6 \im luminosity of the sample 
of Cepheids. 



meaning the condensation radius is significantly large that it is not affected by the pulsation 
and shocks, supporting the predicted results. In all, the dust model discussed implies the 
measured infrared excesses and CSE fluxes can have a very large range, making it more 
challenging to use near-IR interferometry for distance estimates. 



6. Mass Loss and Period Change 



The rate of period change of a Cepheid is a measure of its evolution on the instability 
strip due to the rate of change of effective temperature and luminosity or equivalently the 
change of radius an d luminosity. The period change of Cepheids has been investigated by 
Turner et al.l (120061 ) where the period change is 
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based on the Period -Mean Dens ity relation and the small period dependence of the pulsation 
constant Q oc P 1 / 8 ( Fernielll967l ). This relation assumes the mass of the Cepheid is constant 
with respect to time. Starting with the Period-Mean Density relation, the period change 
can be re-derived to include mass loss, 
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Substitu ting the luminosity and effective temperature for the radius for comparison to the 
result of iTurner et al.l (120061 ). this relation becomes 



P 
P 



4M* 6L* 
'7M, + 7L, 



24 T, 



cff 



7 T, 



off 



(64) 



-37- 



This new relation implies the change of period is due to both evolution and mass loss. Because 
the change of mass of a Cepheid is negative, the mass loss will always act to increase the 
period of pulsation. 



The comparison of the observed and theoretical period change of Cepheids by lTurner et al 



( 120061 ) showed models provide a reasonable fit to the observations overall. There are, how- 
ever, exceptions; one is the large number of Cepheids with a negative period change that 
is smaller in magnitude than that predicted by the models. This difference also appears 
in the comparison of positive period change but the difference is smaller. It is likely the 
models overestimate the rate of evolution, but a lower absolute rate of period change on the 
second crossing of the instability strip could be related to mass loss, which would increase 
the period change. The role of mass loss on period change can be tested by plotting the 
fractional difference {P/P — AM / 7 M)/ (P/P) as a function of the fractional period change. 



6] taken from the literature ( 


Fernie 


1993: 


Turner et al. 


Turner & Berdnikov 


2004 




Turner et al. 


20051) and the 



5|) and the comparison is shown in Figure [TTJ 



Mass loss is a contributing factor to those Cepheids where the fractional difference is differ- 
ent from unity. From Figure [TTJ there are only two Cepheids where the mass-loss rate would 
affect the period change of order 1%, X Cyg and RT Aur; both of which have large mass-loss 
rates and small rates of period change. For the case of RT Aur, it should be noted the rate 
of period change was recently determined to be P = 0.082 s/ yr implying the Cepheid is on 



its third crossing of the instability strip ([Turner et al.l 120071 ). This new result differs from 
the rate of period change used here, P = —0.14 s/yr, but the magnitude of period change 
is still small, and mass loss will still significantly contribute to the rate of period change. 
It can be inferred that mass loss does not play a significant role in changing their periods, 
complementing the result of ITurner et al.l (120061 ). 



The large mass-loss rates correlate to the minimal absolute period change of Cepheids. 
This point is emphasized in Figure [121 where the enhancement of the mass loss by shocks plus 
pulsation is plotted against P/P. This implies mass loss is most enhanced when the Cepheid 
is evolving slowly through the instability strip. This correlation can be produced in several 
ways. For example, when the mass-loss rate is large one might expect the radius of the 
Cepheid to be large, causing the gravity to be lower. A larger radius would also correspond 
to a larger luminosity, and a larger luminosity would mean a lower rate of period change 
according to Equation [6H Another way to get a small period change is for the Cepheid to 
be at the blue edge of the instability strip where both the temperature and luminosity are 
largest for that crossing, meaning the ratios T e g/T e g and are small. At that point in 

the instability strip, the period is smallest, increasing the acceleration of gas due to both 
shocks and pulsation, which could increase the mass-loss rate. The first explanation is seen 
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Fig. 11. — The fractional contribution of mass loss towards the period change for Cepheids 
as a function of period change. Cepheids with a period change not affected by the mass-loss 
rate would fall on the dashed line. Deviations from the dashed line measure how much mass 
loss plays a role. 



in Figure [5] for the Cepheids with large mass-loss rates, but there are no examples of the 
second explanation. The fact that there are few Cepheids from the sample evolving near the 
blue edge may account for this as well as the possibility mass loss may only be extreme very 
early upon entering the instability strip, making it very unlikely to catch one in the process. 

The minimal role of mass loss upon the period change implies the period change is 
primarily due to evolution. Near the blue edge of the instability strip the rate of evolution 
slows as the fractional change of luminosity and temperature is small. Also near the blue 
edge, one might expect the mass-loss rate to be large because the period is lower, causing the 
acceleration due to pulsation and shocks to be larger. For Cepheids on the second crossing, 
the rate of period change is P < 0, but near the blue edge of the crossing, mass loss will tend 
to decrease the absolute value of the rate of period change, and potentially even manage to 
change the sign of the rate of period change. This can be tested by considering a sample of 
Cepheids where the period change is small and consistent with third crossing; the Cepheid 
X Cyg is likely on its third crossing with P = 1.52 s/yr. For X Cyg to be on its second 
crossing its mass-loss rate would need to be M > 2.1 x 10~ 5 M Q /|/r; much too large, in the 
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Name 


P(s/yr) 


Name 


P(s/yr) 


U Aql 


4.29 


Y Oph 


766.82 


r) Aql 


3.24 


BF Oph 


-1.23 


RT Aur 


-0.14 


V350 Sgr 


1.53 


VY Car 


-75.92 


U Sgr 


5.11 


5 Cep 


-0.7 


W Sgr 


3.84 


XCyg 


1.45 


X Sgr 


7.96 


/3 Dor 


12.83 


R TrA 


0.23 


( Gem 


-78.88 


T Vul 


1.05 


W Gem 


-27.66 


SV Vul 


-214.3 


T Mon 


2169.72 







Table 6: The rate of period change for some Cepheids with modeled mass loss. 



range observed for Wolf-Rayet stars. The Cephei d SX Car is another Cepheid consistent 



with being on the third crossing (ITurner et a 



2005J), with P = 0.07 s/yr, period of 4.86 days 



and an inferred upper mass limit of 5.7M (lTurnerlll996l ). This would require of minimum 
mass-loss rate of M = 1.67 x 10~ 6 M & /yr, which is an order of magnitude larger than the 
limit seen in Figure HI It has been shown that it is possible for mass loss to significantly 
affect the rate of period change at certain parts of the instability strip. 



7. Conclusions 

In this study, an analytic method is presented to describe the mass-loss rates of Cepheids 
using as input the global parameters that describe a Cepheid. The derivation is based on 
the method for solving radiatively driven winds with additional energy supplied by shocks 
and pulsation. It is assumed the wind is spherically symmetric, quasi-static and uses the 
Sobolev approximation. The limits of the derivation are explored where in the simplest case 
the mass-loss rate is dependent on the ratio of the pulsation plus shocks and the effective 
gravity. It is shown mass loss can be enhanced by up to three orders of magnitude. 

The analytic method was tested using a set of Cepheids for which the necessary global 
parameters are given in the literature. For reference, the mass-loss rates, based on only 
radiative driving, were also calculated. The energy added from pulsation and shocks in- 
creases the mass-loss rates, in some cases by orders of magnitude, predicting rates in closer 
agreement with those inferred from observations. Furthermore the mass-loss behavior of 
Cepheids as a function of period is described, where the long period Cepheids have mass 
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Fig. 12. — The enhancement of the mass loss, the ratio of the shock/pulsation mass loss and 
the mass loss from radiative driving, as a function of the period change P/P. 



loss dominated by radiative driving but the short period Cepheids have a large range of 
mass-loss rates. The large range of values implies the mass loss for short period Cepheids 
depend on the evolutionary state of the Cepheid. 

The model provides reasonable estimates of mass-loss rates of Cepheids but there are 
uncertainties. The mass is the largest uncertainty, as it is difficult to measure in general, 
so a Period-Mass-Radius relation is used. The value of the radius also has an impact; 
the larger the radius, the lower the gravity and hence the larger the mass- loss rate but 



obser ved radii of Cepheids are precise with an uncertainty of a few percent ( IGieren et al. 



19971 ). There is uncertainty added due to the radiative driving coefficients a and k. The 
mass-loss rate depends on the exponent a, and variations of a can cause the predicted mass- 
loss rates to vary. The value of k is less important because it is small and M oc k x l a . More 
detailed analysis of the mass loss of Cepheids would benefit from calculations of a and k at 
temperatures and gravities consistent with Cepheids over a finer grid. 

For some Cepheids, the mass loss seems to be driven primarily by the shocks generated 
in the atmosphere by pulsation. T he rates cal c ulated are dependent on the accuracy of 



the pulsating atmosphere model by iFokin et al.l (Il996l ). Errors in the model of a Cepheid 



atmosphere with mass M = 5.7M , and T c g = 5750 K will affect other the mass-loss rates 
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will be wrong due to the relative scaling. It would be useful to verify the accuracy of the 
model of the shocks in pulsating atmospheres. Furthermore, models of pulsating atmospheres 
could be used to calculate the mass loss of the Cepheid directly to provide a consistency test 
of the analytic method. 

The shock model may also be tested with observations. Shocks propagating in a stellar 
atmosphere ca use the layers in the st ar to move and this is reflected in the broadening of lines 
in a spectrum. iMathias et al.l (120061 ) analyzed the time varying spectra of X Sgr and found 
evidence for two, possibly three, shocks moving through the atmosphere of the Cepheid in one 
pulsation period. This result suggests spectroscopic observations of a Cepheid with intense 
phase coverage would provide information for understanding shocks in Cepheids. X-ray and 
UV observations are important because shocks moving in the atmosphere will ionize material 



which could be detected as X-ray and UV variability (jEngle et al. 



20061 ). The observations 



would provide a test for any hydrodynamic model of shocks in a Cepheid atmosphere. 

It would be interesting to apply the theory of pulsation and shock induced mass loss 
to first overtone Cepheids and other radial pulsating variables. As first overtone pulsation 
occurs at an earlier stage of pulsation than fundamental pulsation, analysis of mass loss in 
first overtone pulsating Cepheids would shed more light on the contribution of mass loss on 
the problem of the Cepheid mass discrepancy as well as contribute t o the understanding 
of Cepheid evolution. However, the shock model of iFokin et al.l (119961 ). that underlies this 
analysis, is valid only for fundamental pulsation and thus cannot be realistically applied to 
first overtone pulsation. With an appropriate shock model, shock induced mass loss can be 
investigated with this theory. 

The dust model proposed here provides a consistent match with existing observations. 
In particular, the calculation of the ratio of the l umin osity at wavelengths 25 fim and 12 \im 
exhibits the same behavior as shown by iDeasyl (119881 ). The results, however, are dependent 
on the choice of dust parameters such as mean grain size. New observations from Spitzer 
and at sub-millimeter wavelengths would provide constraints on the properties of the dust 
in the wind and on the mass-loss rates of Cepheids. Sub-mm observations would also 
provide information on the age of Cepheids. Consider RS Pup with four observed shells, 
each associated with a particular crossing on the instability strip; a very strong sub-mm 
excess may indicate a large amount of cool dust in these shells. Thus the amount of sub-mm 
excess may trace the mass-loss history of the Cepheid and its age. 

The connection between mass loss and evolution has been explored for the observed 
Cepheids. It was shown the period change of a Cepheid is dependent on the mass-loss 
rate at that time. However, for all the Cepheids with period changes quoted, the mass-loss 
rate provides an insignificant contribution, and, furthermore, the Cepheids with the largest 
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calculated mass-loss rates have the slowest absolute value of period change. This connection 
needs to be explored in greater detail with a larger set of observations. 

Mass loss in the instability strip is an important concept to understand if one is to 
determine the physical structure and the cause of the mass discrepancy of Cepheids. It is 
also important to constrain mass loss for observations in the near-IR and longer wavelengths 
if the period-luminosity relation is to be precise to a few percent, which is a motivation of 
interferometric observations. The analytical model of mass loss of the Cepheids investigated 
here imply shocks and pulsation play a strong role driving the wind and increase the mass- 
loss rate by orders of magnitude. 

HN would like to thank John Percy for careful reading and comments, and Tom Bolton 
for teaching the course on stellar winds that inspired this work. The authors would also 
like to than the anonymous referee for helpful comments. The research was funded by the 
Walter John Helm OGSST and the Walter C. Sumner Memorial Fellowship. 
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